Resonant tunnelling and quenching of tunnel splitting in Wess-Zumino nanospin 

systems 
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We investigate the energy spectrum of the biaxial spin systems with magnetic field along the hard 
anisotropy axis by using the complex periodic orbit theory. All important features of the system 
appearing in whole energy range, such as oscillations of level splittings due to Wess-Zumino effect 
and their absence at higher magnetic field etc., can be completely understood within this semiclassical 
scheme. We find out that the fields at which the tunnelling quenches do not shift at higher energy 
levels and the absences of the quenching at higher magnetic field have their origin in an exact 
coincidence of the quenching field with the field of resonant tunnelling. Based on the result, we 
propose that the complete cancellation of quenching with resonant tunnelling would be a general 
property of Wess-Zumino tunnelling systems. 

PACS numbers : 75.45.+j, 03.65.Bz, 03.65.Sq 
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Recently, the oscillations of the level splitting of the 
nanospin systems with magnetic field along the hp 
axis has been studied by many authors theoretically 
and experimentallyu. This oscillation phenomenon is 
semiclassically understood as a result of the interference 
between instanton and anti-instanton in the Euclidean 
spaceo. However, although this instanton method gives 
exact evaluation of the oscillation of the ground level 
splitting, it does not give a satisfactory explanation of 
how the quenching disappears at higher magnetic field 
and how higher level splitting behaves with magnetic 
field. 

On the other hand, there have been many works 
about complex periodic orbit theory for the tunnelling 
systemsBia This theory is basecLon the extension of well- 
known periodic orbit theoryBliil by including tunneling 
paths going through energy barrier. The tunnelling path 
have a minimal Euclidean action with, a given eneijgy, so 
they can be regarded as a 'bounceO, 'instantonl!3, or 
'periodic instantonO. The advantage of complex peri- 
odic orbit theory is that for tunnelling systems it gives 
energy levels and their splittings throughout whole en- 
ergy range. It is worth while to see how the geometrical 
phase can be incorperated with the complex periodic or- 
bit theory, because, to our knowledge, for the systems 
with geometrical phase the theory has not been applied 
so far. 

The purpose of this Letter is to apply the complex 
periodic orbit theory to the nanospin system showing ge- 
ometrical phase effect, and to give semiclassical explana- 
tion for the existence of the level-splitting oscillations at 
low magnetic field and the absence of the oscillations at 
higher magnetic field for excited states as well as ground 
state. As will be shown clearly, we find out that the 
mechanism of the disappearance of quenching is the exact 
cancellation of quenching with the divergence of resonant 
tunnelling, which is very novel feature not reported so far. 



This strong correlation between quenching and resonant 
tunnelling seems to be inherent property in Wess-Zumino 
tunnelling systems. 

Let us start with the Hamiltonian of system with mag- 
netic field along the hard axis, 



H = KxSi 



g^ B HS z 



(1) 



where K\ and Ki are anisotropy coefficients with K\ > 
K%, and g and fj,B are the gyromagnetic ratio and the 
Bohr magneton, respectively. The effctive Euclidean ac- 
tion for this Hamiltonian is obtained by expressing it in 
(0, 4>)- representation and integrating over (k3, 
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where Sw denotes the Wess-Zumino term which is ex- 
pressed as 



Sw = S / dr4>(l — 



1 - A sin^ 



-)• 



(3) 



and the position dependent mass and effective potential 
are 



m(<p) 



1 



1 - A sin^ 



n^) = ^sin 2 0(l- f, ). 
I 1 — A sm m 



(4) 
(5) 



Here we put A 



#a and h 



Note 
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that the geometrical phase factor Sw does not change 

under analytic continuation of Euclidean time into real 

time, i.e., r — > it. This means that one should consider 

the phase effect of Sw for the path in real space in the 

same way as for the tunnelling path in the Euclidean 

space. This— point of view is somewhat different from 

that of Ref.Ej where only the phase effect of tunnelling 
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path is taken. The lack of the contribution of real path 
to the phase effect has created unreasonable oscillation 
of level splitting with temperature. 

Now, we apply complex periodic orbit theory to this 
effective Hamiltonian system with Wess-Zumino term. In 
the regime of h < hi = 1 — A, as depicted in Fig. 1 (a), 
the potential shows a simple barrier between <fi = and 
7r. Let us first consider the energy range of E < Eb, 
Ei, being the energy of the barrier top. A periodic orbit 
starting from cf> = <f>\ can be considered as a combination 
of four elements; (tunnelling segment + tunnelling seg- 
ment ), (real segment + real segment), (tunnelling seg- 
ment + real segment), and (real segment + tunnelling 
segment). With a prescription for Maslov index for tun- 
nelling segment in Ref.cl, we can write semiclassical trace 
of Green's function in terms of a periodic orbit sum as 



T 



I 

+ (e 
+ (e 



m=0 N=0 
— VgtW— i 

iW-i%\f„-6+i 



r)(e" 



<-' + j#\2 _|_ / iW-if 

■ya e im-K 

ia irmr~\ N 



■ \ —VOL K 

)e e 



(6) 



where T is the period of the real orbit in the well, W is the 
action for the real segment, i.e., W(E) = f_\ d<j)m(<j))^, 
9 is the action for the tunnelling segment, i.e., 9(E) — 
* 1 d(f>m(<j>)^!p, and a — Sw(l — ) is the geomet- 
rical phase for an element corresponding to tt increment 
in <f> axis. The first two terms correspond to the ele- 
ments (tunnelling segment + tunnelling segment) and 
(real segment + real segment), respectively. These ele- 
ments represent oscillatory motions confined in the wells 
and barriers, which do not enclose z-axis and thus have 
no geometrical phase factor. While, the next two terms, 
(tunclling segment + real segment) and (real segment + 
tunnelling segment) represent half rotations in the direc- 
tion of increasing (j> and decreasing <f>, respectively, so 
that they have the geometrical phase factor e ta or e~ la 
according to the rotation direction. The contributions 
of nonperiodic orbits ending at 4>i + tt in the combina- 
tions are removed by introducing m sum. Using the same 
approximation as Refill, we can write the sum as 



9(E) 



^-f E - (Ek ± y COS ( 



(J) 



where Ek is the EBK(Einstein-Brillouin-Kcller) energies 
determined by 2W(E k ) = 2n(k + 1/2), k being integer, 

and T — ^j|p. Since the poles of g(E) indicate the 
eigenvalues of the system, we can see the oscillation of the 
energy splitting and the quenching fields are determined 
by a(h) — ir(n + 1/2), n being integer. It is noted that 
these quenching fields do not shift in excited energy level 
splittings. 

Now, consider the regime of h > hi. In this case the 
barrier top becomes local mininum and then equivalent 
two barriers and new meta-stable well appear as shown 
in Fig. 1 (b). Eq.(Q) is still valid in the energy range of 
E < E m . Consider the energy range of E m < E < Eb, 



E m being the local minimum energy. In this structure 
one can expect that the resonant tunnelling may happen 
due to the states in the meta-stable well. As will be 
shown below, however, the resonant tunnelling does not 
take place, since it compensates the quenching of the 
level splittings. As a result, there is no more oscillations 
of level splittings. The periodic orbit sum can be done 
simply if one attaches a small correction 



1 - 



(8) 



to each segment and replace the contribution of tun- 
nelling segment by 



-26 m +iTT p iW m -in/2 



1 



,2iW m -i 
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where W m (E) - ^ d^m^)^ and 9 m (E) = 

d(jym(4>)^. The small correction, Eq.(||), comes from 
the consideration of possible paths, leaving from <pi or 
4>i + 7r and returning back after some oscillatory motion 
in the meta-stable well, and their repetition between two 
segments. In fact, this small correction does not play any 
important role under the approximation used in getting 
the final g(E). The contribution of the new tunnelling 
segment, Eq.(^|), shows all possible oscillatory motions in 
the meta-stable well between two equivalent tunnelling 
subsegments. Then, the trace of Green's function can be 
written approximately as 
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Here, it is noted that the geometrical phase effect cos a is 
modulated by the action of meta-stable state. This effect 
is the main point in this Letter and will be discussed in 
detail soon. Note that in the limit E — > E m , Eq.(|10|) 
gives smaller energy splittings than those given by Eq.(M, 
i.e., shows a discrepancy of factor 2. This fact leads to a 
discontinuity at E = E m in energy splittings (see Fig. 3). 
To remove this discrepaHjCT at E — E m , one should use an 
uniform approximationQO which may also give a smooth 
behavior near Eb as well as near E m . However, we do 
not use the uniform approximation because the simple 
complex periodic orbit theory used here is sufficient for 
our purpose of this Letter. 

If the energy is larger than Eb, the particle motion 
is simple rotation and the poles of the trace of Green's 
function appear when 



Srot(E) ± 2a = 2irk, k — integer, 



(11) 



where the action S ro t(E) = d<pm(<f>)^ and this con- 
dition gives EBK energies Ek for E > Eb- In this energy 
range the energy splitting does not come from tunnelling, 
but rather from the Zceman effect. Indeed, this simple 
semicalssical evalution of EBK energies gives exact eigen- 
values for uniaxial spin system of H = S% — HS Z which 
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is not a tunnelling system but an example of simple rota- 
tional orbit. In this uniaxial spin system the mass is 1/2, 
the potential is —H 2 /4, and a = ir(S — H/2), and then 
from Eq. ( |TT| ) the EBK energies is given as Ek = k 2 — kH 
where k is integer or half integer according to the spin- 
value S. 

All actions, W(E), 6(E), W m (E), 9 m (E), and S rot (E) 
are analytically calculated by the complete elliptic inte- 
grals of the first and third kindlll (see Table 1). 

In Fig. 2 (a) the whole energy spectrum with h are 
shown for S = 10 case. The solid line and the dashed 
line denote the barrier energy Eb and the local minimum 
energy E m , respectively, with h. In order to avoid com- 
plicated view, we draw only EBK energies for E < Eb 
without showing the splitting of the EBK energies. For 
the comparison we draw the result of direct diagonaliza- 
tion for the same system. In this direct numerical cal- 
culation we take the energy such that E = at the well 
minimum. In fact, the effective potential in Eq.(^|) has 
been derived after the same energy translation. We can 
see that the entire behaviour of energy spectrum with 
h is very similar to the result of direct diagonalization 
except for the region of E ~ Eb- As mentioned before, 
an uniform approximation would give a consistent ex- 
planation for this transition region. In our semiclassical 
evaluation of the energy spectrum the physics in the en- 
ergy range of E > Eb is very different from that in the 
range of E < Eb] the former is the Zeeman splitting and 
the latter is tunnelling splitting. However, in quantum 
mechanical result these seem to be mixed around E ~ Eb 
and no abrupt change at E — Eb is shown. 

Fig. 3 shows the lowest three level splittings, AE Q , 
AEi, and AE%. Here, the solid lines describe the results 
of complex periodic orbit theory(A£"|), while the dashed 
lines show the result of direct diagonalization(A£^). It 
should be emphasized that in the evaluation of the ac- 
tions related to tunnelling segment, i.e., 0(E), W m (E), 
and 6 m (E) , the energy value used is not the EBK energies 
E k , but El determined from the condition 2W(El) = 
2ixk. This choice of energy is consistent with usual in- 
stanton method where the ground state splitting is rep- 
resented in terms of Euclidean action at the well mini- 
mum energy E m i n (action of instanton) instead of that at 
the ground state energy Eq. The discrepancies between 
the quantal results (AE^) and the semiclassical results 
(AE-) would be resolved provided that the quantum fluc- 
tuations around the tunnelling path are taken into ac- 
count. It is very interesting to see that the quenching of 
the level splitting does not exist at higher h, and for AEq 
the number of quenching is S and AEi has (S—i) quench- 
ings in the range of < h < 1. This can be understood 
if one knows the role of cosW m in Eq. (|To|) . The factor 
1/ cos Wm diverges at 2W m = 2n(n+ 1/2), n being inte- 
ger, showing the resonant tunnelling, which means that 
if El coincides with the EBK energy E™ eta in the meta- 
stable well the level splitting becomes divergent. The 
fact that in the range of E m < Ei < Eb the splittings 
AEi (i — 0,1,2,3) do not show any quenching and any 
resonant tunnelling implies that the quenching fields are 



same with those of resonant tunnelling. This exact coin- 
cidences are shown in Fig. 4. The dashed lines denote the 
EBK energies in the meta-stable well, ij™ eta . Note that 
the cross points with Eq(= 0), E{, and E% (the fields 
of resonant tunnelling) are located exactly on quenching 
fields given by the condition cosa(/i) = 0, which is the 
mechanism for the disappearance of the level-splitting 
oscillation. 

In conclusion we examine the energy spectrum of the 
quantum tunnelling Wess-Zumino spin system by means 
of the complex periodic orbit theory. All important fea- 
tures in the spectrum can be understood within this semi- 
classical scheme. Especially, we find out that the absence 
of the oscillation of the level splittings at higher magnetic 
field has its physical origin in the exact cancellation of 
quenching with the resonant tunnelling phenomena. It is 
more reasonable to conclude that this cancellation is an 
inherent property of Wess-Zumino tunnelling system. 

We thank Stephen Creagh and C.-H. Kim for useful 
discussions. 
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FIG. 1. The effective potential Eq. @. ( a) h < hi case 
(A = 0.7 and h = 0.1). The solid and dashed arrows de- 
note the real segments with action W and the tunnelling seg- 
ments with Euclidean action 9, respectively, (b) h > hi case 
( A = 0.7 and h = 0.5). The solid arrows in the meta-stable 
well represent the path with action W m , while the dashed 
arrows do the tunnelling subsegments with Euclidean action 

dm- 



FIG. 2. The whole energy spectra for the Hamiltonian Eq. 
(0) with spin value 5" = 10 and A = 0.7. (a) The result of 
the present semiclassical theory. The solid and dashed lines 
represent the energy of barrier top Eh and the energy of the 
local minimum E m , respectively. For the energy range of 
E < Ei we draw only the EBK energies, (b) The result of 
direct diagonalization for the Hamiltonian Eq. (Fh. 



FIG. 3. The energy level splittings AE , AEi, and AE 2 
in the case of S = 10 and A = 0.7. The solid lines denote the 
result of the complex periodic orbit theory (A-Ef) and the 
dashed lines do the result of direct diagonalization (AEi). 



FIG. 4. The EBK energies in the meta-stable well 
E™^ 11 (dashed lines) and E£ (solid lines). The cross points 
imply the resonant tunnelling. The cosine function is 
0.06 cos a(h). Note that the quenching fields given by 
cos a(h) = are located on the same h with the fields of 
resonant tunnelling. 



TABLE I. Analytic expressions for various actions. K(k) and il(a 2 ,fc) are the complete elliptic integrals of first and third 
kind, respectively. 



W(E) = 
0(E) = 
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[A(B. - i)II(a?, ki) + (1 - S-)n(ai ki) 



A — X 



K(ki)] 



W m (E) = 
9 m (E) 



S /X(l-B_)B + L * 

1 r 2XE-H 



W ^ T) n(a§, hi) + B-B(al k 2 ) - ^K{k 2 )] 
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X /X(l-B_)B + A 

Srot{E) = iA [n(a|,fe 5 ) - (1 - XB-)U(al ,ks)] (in the case of real B- 
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